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Abstract: We obtain explicitly the renormalization group equations for the quark mass 
matrices in terms of a set of rephasing invariant parameters. For a range of assumed 
high energy values for the mass ratios and mixing parameters, they are found to evolve 
rapidly and develop hierarchies as the energy scale decreases. To achieve the experimentally 
observed high degree of hierarchy, however, the introduction of new models with specific 
properties becomes necessary. 
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1. Introduction 

A long-standing problem in particle physics is to have some theoretical insights regarding 
the plethora of parameters contained in the mass matrices of quarks and leptons. These ma- 
trices arise from the coupling of the Higgs bosons to the fermions, and are not constrained 
by any known principle. The coupling gives rise to ten parameters each in either the quark 
or the lepton sector. These seemingly arbitrary quantities do have some regularities. In the 
quark sector, both the mass ratios and mixing parameters exhibit rather large hierarchies. 
On the other hand, while the lepton masses are hierarchical, the neutrino masses tend to 
be more degenerate with some large mixing angles. When one tries to sort out possible 
clues contained therein, we need also keep in mind that these parameters are all measured 
at low energies, and it is necessary to bring renormalization effects into the picture. As 
the energy scale changes, one expects the pattern of regularity to evolve according to the 
renormalization group equations(RGE). Thus, we might entertain the hope that a simpler 
picture will emerge at high energies. The analysis of the RGE of the quark mass matrices 
has a rich literature and a long history [1-9]. However, general conclusions are not easily 
available. For one thing, the RGE are simple when formulated in full mass matrices. This 
means that there are a large number of superfluous degrees of freedom which must be 
stripped away to get at the physical variables. However, the choice of these variables are 
not unique. In the literature, the extant RGE turn out to be very complicated and highly 
nonlinear. As a result, while their low energy behavior is generally known, very little can 
be said of their intermediate and/or high energy behaviors. 
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In this paper we propose to write the RGE in terms of a set of rephasing invariant 
mixing parameters which were introduced recently. The resulting equations are simpler 
than those given in terms of other parameters and are amenable to a general analysis. At 
low energies, it is found that the RGE are close to a fixed point when we put in the physical 
values observed experimentally. This result is well-known in the literature, so there is very 
little evolution of the parameters in the low energy region. If one were to extrapolate this 
low energy behavior to all energies, then one might conclude that renormalization effects 
are altogether unimportant. However, this extrapolation is based on the assumption that 
the RGE are exact for all energies and that the initial values are precisely known. A more 
likely scenario is that what we have at low energies is an effective theory which comes 
from some new theory valid at high energies. It is therefore more appropriate to evolve 
the RGE from high to low energies. To do this we will assume that the parameters are 
generic, instead of being hierarchical, at high energies. It is seen that the set of RGE 
does give rapid running for a range of initial values. This suggests an "infrared hierarchy" 
scenario, viz., low energy hierarchy comes from renormalization evolution for a range of 
generic values at high energies. To account for the very large hierarchy (up to e~^^ ~ A^, 
A ~ 0.23) observed experimentally, however, it is necessary to invoke new theories different 
from the standard model(SM) or its minimal supersymmetric extension(MSSM). In the 
absence of concrete models we mimic their effects by changing some parameter in the 
known RGE. It is found that, for appropriate choices, all mixing parameters evolve very 
rapidly and infrared hierarchy is a viable scenario. However, mass ratios do not develop 
a strong enough hierarchy to match the observed values. Obviously, more work needs to 
be done and it is necessary to construct explicit models to realize these ideas. We hope to 
report on our progress in a future publication. 

This paper is organized as follows. In Secj2|, we review and summarize the (x, y) 
parametrization introduced earlier. The RGE for {x, y) and mass ratios are given in Sec. |3|. 
Section-]^ is devoted to the exploration of solutions for these equations which are compatible 
with the idea of infrared hierarchy. After some concluding remarks in Sec. ^, we discuss 
in Appendix ^ detailed properties of the matrices which enter the RGE. Finally, brief 
summaries on two-loop renormalization are included in Appendix ^ 

2. Rephasing Invariant Parametrization 

To establish our notation and for completeness, we now review briefly the rephasing in- 
variant parametrization introduced earlier [10, 11]. 

Three flavor mixing is described in terms of a 3 x 3 unitary matrix Vij, where i,j = 
1, 2, 3. Without loss of generality, we may impose the condition det V = +1. Then, there 
are six rephasing invariant combinations 



where {i,j,k) = cyclic permutation of (1,2,3) and J is the Jarlskog invariant [12], which 
turns out to be the common imaginary part of all Tjjfc. We define 



^ijk — yiiV2jV^k — Rijk 



— iJ, 




(Xl, X2, X3; yi, 2/2, ys) — (-Rl23) -^231) -R312; -R132, -R213, -R321) 
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It is found that they satisfy two constraints 



det V = {xi+X2 + X3) - {yi + 2/2 + 2/3) = 1, 



(2.3) 



X1X2 + X2X3 + xzxi = yiy2 + 2/22/3 + 2/32/1, 



(2.4) 



leaving four independent parameters to describe the mixing. In addition, there is a simple 
relation 



J = X1X2X3 - 2/12/22/3- 
These parameters are all bounded by ±1, i.e., 



-1 < (a;i,2/j) < 1, yj<Xi. 



(2.5) 



(2.6) 



They are related to \Vij\'^ by 



W = 



|Vii|2 \Vi2\^ |Vi3|2 
1^21^ IV22P 1^23!^ 
|^31P |^32P \Vss\^ 

Also, the matrix of the cofactors of W, with w'^W = (det W)I, is given by 



Xl -yi X2- 2/2 3^3 - 2/3 

iC3 - 2/2 ici - 2/3 a;2 - 2/1 
y a;2 - 2/3 - 2/1 a;i - 2/2 



(2.7) 



It; 



Xl + 2/1 X2+ 2/2 X3 + 2/3 
X3 + y2 Xi+ 2/3 X2 + 2/1 
y X2 + 2/3 3^3 + 2/1 3^1+2/2 



(2.8) 



The matrix w appears repeatedly in the RGE which will be presented in the following 
section. Experimentally, the elements of the quark mixing matrix Vckm are well measured 
and they exhibit striking hierarchies: 



(xi,X2,X3;-2/i,-2/2,|y3|) = (1,A^A^+A^+A^A«). 



(2.9) 



It should be emphasized that, unlike the usual CKM parametrizations, which contain terms 
of order A, here the hierarchy is in powers of A^ only, which is also the hierarchy pattern 
in quark mass ratios. We add that, since we can choose any set of four of the {x,y) 
variables to parametrize Vckm-, including 2/3 (which is very small) in such a set implies 
that the physical Vckm is effectively described by three parameters only. Thus, there 
are correlations amongst the familiar four parameter sets in Vckm- These approximate 
relations [11] are 

Sl3 = C5S12S23, (2.10) 



\Vus?\Vcb? = \Vtd? + \Vuh?. 



All of these are well satisfied by existing data [13]. 
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c 


b 


a« 




SM 


-3/2 


3 


2{T - Gu) 


2(r - Gd) 


MSSM 


+1 


6 


2(3TriV4 - Gi) 


2(3TrMrf + TrMe - G^) 


THM 


+1/2 


3 


2(3TrM„ - 


2(3TrMrf + TVMe - G^) 



Table 1: Coefficients in Eqs.(|]|, U). 



3. Renormalization Group Equations 

The RGE for quark mass matrices have been obtained and studied for a long time. Al- 
though they are simple when expressed in matrix form, they become complicated and 
highly nonlinear if one writes them in terms of the parameters commonly used for the 
CKM matrix. In fact, their complexity has been the major obstacle preventing one from 
drawing general conclusions so that most analyses are confined to low energies and their 
extrapolations [6-9]. 

In this paper we will reformulate the problem using the (x, y) parameters summa- 
rized in the previous section. The resulting equations turn out to be simpler and not so 
formidable. In addition, when we write down the evolution equations for some simple func- 
tions of the mass ratios, they take very similar forms to those of the {x, y) parameters. This 
bolsters the idea that the observed hierarchies in the mass ratios and the (x, y) parameters, 
which so closely resemble each other, may both be related to renormalization. 

We begin by citing the one loop RGE for the mass (squared) matrix of the u-type 
quarks, = YvYu, and that of the d-type quarks, = Y^Y^, where Y denotes the 
Yukawa coupling matrices of the Higgs boson to the quarks [5-7]. 

= auMu + bMl + c{M„, MJ, (3.1) 

^Md = adMd + bMl + c{M„, Md}. (3.2) 

Here, Si = lOvr^^ and t = \n{ii/Mw), where H is an energy scale and Mw is the W boson 
mass. The values (a^, ad,b, c) are model-dependent and we list them in Table |l|, for the 
SM and the MSSM with tan/3 = 1, as well as the two Higgs model(THM), where one 
Higgs couples to u-type quarks, and the other to d-type quarks and the leptons. Here, the 
notations used in Table || follow, 

r — 17^2 I 9^2 I Q„2 r' — 1^2, 9_2 |o2 

Gi = M + Sgi + fgl G^ = + 3<7i + fgl (3.3) 
T = TV(3M„ + 3Mrf + Me), 

where gi, i = 1,2,3, are the usual gauge coupling constants. 

For our purposes it is convenient to start with the explicit evolution equations of 
the eigenvalues of the mass matrices and those of the CKM matrix elements [6, 7], using 



Eqs.(3.1, p. 2D . For the eigenvalues, we have 

= f!K + bf! + 2c^h]\V,,\% (3.4) 
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and 



= h] [ad + hh] + 2c ff\Vi, 



(3.5) 



where ff and /i| denote the eigenvalues of and M^, respectively. The explicit equations 
for the CKM matrix elements are 



where we have defined 



Fikh^VieYkuVkj + y, HjkfmVmk^mjVik 



Fik 



m,k^j 



n-fi' 



H. 



jk 



h] + hi 



(3.6) 



(3.7) 



Eq.(|3.6|), as it stands, is not rephasing invariant. To it one could append terms with 
i = k (or j = k) on the right hand side [6, 14]. They come from (purely imaginary) 
diagonal elements of {UdW /dt), where U diagonalizes Mu{ot M^). A rephasing transfor- 
mation, U — > (exp ia(t))C/, where a{t) is a t-dependent, diagonal, phase matrix, yields 
{UdW/dt) e'°'{UdUydt)e-'''' - ida/dt. Thus, the diagonal elements of (UdUydt) are 
rephasing dependent, and Eq.( |3.6|) should only be used to compute the evolution of rephas- 
ing invariant combinations of Vij, such as \Vij\'^ or {xi,yj). 

Using these equations, while keeping in mind that the condition det^ = +1 implies 
relations such as V-^i = V22V33 — V23V32, etc., we can work out the RGE for the mass ratios 
and the mixing parameters in terms of the (x, y) parameters. We find, after some algebra. 



-^Xi/C= (A/23, A/31, A/12) 




(A/123, A/131, A/112) Si 



/ i?23 \ 



-%,/c= (A/23, A/31, A/i2)4 



'23 

where we have defined the differences of mass squared 



+ (A/l23,A/l3i, A/112 )5,' 




Af,j = f^-ff, Ah, 



hf - hi 



(3.8) 



(3.9) 



(3.10) 



The matrices Ai, Bi, A'- and B'- are given in Table |2[ Note that all of the elements of 
these matrices are bounded [0 < (IxjXjl, \xiyj\, \yiyj\) < 1]. As will be discussed in detail 
in Appendix these matrices satisfy consistency relations which guarantee the identities 



^'^ixi - yi) = 0, ^^{xiXj -yiyj) 



i<j 



Also, by calculating ^(^1X2X3) and ^(2/12/22/3) (see Eqs.( A.6 
evolution equation for J^: 



0. (3.11) 
1)), we obtain a simple 



-2cj2(A/^ -wH + Ah^ 



w 



(3.12) 
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where we have used a compact notation for the column matrices: 
{Af,Ah,F,H} 





( A/23 \ 


( A/123 \ 


( F2A 


( ^''^ ] 


< 


A/31 , 


A/131 , 








\Ml2) 


A/112 / 







(3.13) 



Eq.( |3.12 ) agrees with previous results. 
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Ai 



Xi 



^ yi X2 X3 
X3 2/3 X2 
\ X2 X3 y2 



+ 



^ yixi yzy2 ^22/3 
Viy2 ysxi y2yi 
\yiy-i ysyi y2Xi 



Br 



Xl 



^ yi x^ X2 
X2 ys X3 

\ X3 X2 2/2 



+ 



yixi y3y2 2/22/3 
2/12/2 y-ixi 2/22/1 

\ 2/12/3 2/32/1 y2Xi 



X2 



^ Xl 2/2 3^3^ 




2^3 Xl 2/1 




\ 2/3 2^3 Xl j 





2/32/1 y2X2 2/12/3 

2/32/2 2/22/3 2/12:2 







^ Xl X3 2/3 ^ 




) 


X2 


2/2 2;i X3 








\X3 2/1 a;i y 





2/22/1 2/12/2 2/3a;2 

2/22;2 2/12/3 2/32/1 



a^3 





2/3 ^ 




2/2 Xl 


X2 


K 


\X2 2/1 


Xl j 





2/22/1 yi2/2 2/3x3 

2/2X3 2/12/3 2/32/1 



X3 



^ Xl 2/2 X2 
X2 Xl 2/1 
\ 2/3 X2 Xl 



+ 



^ 2/32/1 2/2X3 2/12/3 
2/32/2 2/22/3 2/ia;3 
\ 2/3X3 2/22/1 2/12/2 



2/1 



^ Xl 2/2 2/3 ^ 




2/2 2/3 X2 




\ 2/3 X3 2/2 y 





X12/1 X3X2 X2X3 
X1X3 X3X1 X22/1 







Xl 


2/2 


2/3 ^ 




) 


2/1 


2/2 


2/3 


X3 








\2/3 


X2 


2/2 y 





X12/1 X2X3 X3X2 
X1X2 X2X1 X32/1 



2/2 



'^2/1 a;2 


2/3 ^ 




X3 2/3 


2/1 




\2/3 2/1 


Xl y 





X3X1 X22/2 X1X3 
X32/2 X2X1 X1X2 



2/2 



2/1 2:3 2/3 
X2 2/3 2/1 
\2/3 2/1 Xl 



+ 



( X2X1 X32/2 X1X2 
X22/2 X3X1 X1X3 
\x2X3 X3X2 X12/2 



2/3 



'^2/1 2/2 X3 
2/2 Xl 2/1 
\ X2 2/1 2/2 



+ 



'^X2Xi X1X2 X32/3 
X2X3 X12/3 X3X2 

\ X22/3 2:1X3 X3X1 







^ y\ 2/2 X2^ 




) 


2/3 


2/2 Xl 2/1 








i^X3 2/1 2/2 y 





X3X1 X1X3 X22/3 

X3X2 X12/3 X2X3 



Table 2: The matrices Ai, Bi, A[ and B[ used in Eqs.(3.8, 3.9) 



-7- 



We now turn to the RGE of mass ratios. Prom Eqs.(3^, |3.5| ), it is straightforward to 
obtain 



( lni?23^ 

lni?3i 
\lni?i2 1 


1- 


/A/23^ 
A/31 

VA/i2y 




' In ^ 

In r3i 
\lnri2 / 




( A/123 ^ 

A/131 
V A/112 ) 


■{ 



+ y2 xi+ ys X2 + yi 



X2 + y2 Xi+ 2/3 X3 + yi 




(3.14) 



(3.15) 



Here, we have defined ratios of masses as 



Ri 



fi/fj 1 



hi/hj. 



We summarize this by writing Eqs.(3.14,3.15) in the form 

^lnR= ^Af + cwAh, 



&lnr = -Ah + cw^ A/, 
where we used Eq.( 3.13| ) and the definition 



(3.16) 

(3.17) 
(3.18) 







/lni?23\ 


^ In r23 \ 1 


(Ini?, Inr) = < 




lni?3i , 


lnr3i . 






\lnR12 J 


i^lnri2 J J 



(3.19) 



In going from eigenvalues (/j, hi) to their ratios, we find a set of much simphfied RGE 
which only depends on the mass differences and the (x, y) parameters in the combination 
w and w"^ , defined in Eg. ( p. 81) , a fact which seems very interesting but not understood. 
The RGE for the mass ratios, Eqs.(|l|, |1|), can be brought into a form in conformity 
with those for the {x,y) parameters by considering sinh(lni?jj) = (/? — fj)/2fifj, with 



coth(ln i?jj). Note that in the hierarchical limit sinh(lni?j 



for /i »/j. 



Let us define 



Qij = ln[sinh(lni?jj)], qij = ln[sinh(ln rj^)]. 



(3.20) 



We find 



i^(Q23,Q31,Ql2) 



|(A/23, A/31, A/12) 



/ F23 
F31 
\ F12 



+c(A/l23, A/131, A/112 )w^ 



/F23 \ 

F31 
V F12 / 



(3.21) 
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^(923,931,912) 



K A/123, A/131, A/l 



12 



/ if23 \ 
^31 

V Hi2 J 



I Hoz 



(3.22) 



+c(A/23,A/3l,A/i2)ti; 



V 



'23 

Hu 
Hi2 



The evolutions of Qij and qij are thus similar in form to those of {xi,yj) and J^. In 
fact, using these equations we readily recover the elegant result for the evolution of the 
CP— violation measure [15], which is a product of J and the quark mass-squared differences. 



^in[jn(A/i,-) n(^^fc^)/(n/^n^')] = ^Tr(M„ + 



(3.23) 



where b is given in Table |l[ Eq.( p.23|) agrees with Ref. [15] when we use the RGE of det 
and det M^, given by 



^(IndetM^) = 3a„ + 6TrM„ + 2cTrMd, 



(3.24) 



^(In det Md) = 3arf + bTvMa + 2cTrM^,. (3.25) 

The ten independent equations, which are contained in Eqs.(U, |]2l|, |3]2|, |]2|, 
form a complete set of RGE for the ten physical parameters in and M^. Although it 
does not seem feasible to analytically[16] disentangle these coupled, non-linear, differential 
equations, their dependence on the relevant variables are rather simple. So their analyses 
are certainly less demanding compared to the corresponding equations in the literature, 
which were written in terms of other parametrizations. 

We offer some general remarks. Let us first observe that the RGE for the (x, y) 
variables, Eqs.(3.8, p.Op , have a fixed point at the parameter set (or permutation thereof): 
xi = 1, Xi = yj = 0, i 7^ 1. The physical values for {xi,yj) are thus very close to this 
fixed point, and little evolution is expected of these variables within the low energy region. 
As for the evolution of mass ratios, note first that for hierarchical ratios, Eqs.( |3.14 , 3.15| ) 



coincide with Eqs.( p.21 , 3.22| ), with their fixed point at A/^ = Ahij = 0. Thus, in the 



physical low energy region, with w = w —I, and the only appreciable physical value 
being ~ 1, again the mass ratios do not change much. (Numerically, the evolution over 
a region near t = 0, say At ~ 10, !^Rij ■ At/Rij would be no more than 10%.) We can 
safely conclude that, in the low energy region, renormalization effects are small and little 
deviations are expected of the measured (x, y) values and mass ratios. This conclusion, 
that the low energy physical values are close to a fixed point of the RGE, confirms the 
known numerical analyses given in the literature. One could extrapolate these results to 
high energies. But to do so accurately assumes that not only the RGE are exact, but also 
the parameters are very precisely known. We would rather examine these equations from 
the high energy end and study their behavior which might impact low energy physics. The 
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structure of the RGE suggests that one can have rapid evolution if 1) masses are large; 2) 
the {x,y) values are large; 3) degeneracy so that Fij and/or Hij are large. Under theses 
initial conditions it is possible to realize the "infrared hierarchy" scenario, namely, the low 
energy hierarchies originated from RGE which quickly drive the relevant quantities toward 
the observed values. As the energy scale decreases, the high energy theory is taken over by 
the low energy models. However, they will have little effect on the established hierarchies. 



4. Numerical Results 



In the absence of analytical solutions for the RGE, we turn to a numerical analysis of these 
equations. As we argued in the previous section, while the low energy physical values are 
close to those of a fixed point of the RGE, it is not clear how fast they are approached from 
presumed, more generic, high energy values. The structure of RGE suggests that rapid 
evolution is possible if all of the {xi,yj) parameters are of order one and that the masses are 
large and nearly degenerate. In addition, the constants b and c, which are characteristic of 
the underlying model (Table [l|), are important factors in determining the evolution speed. 

Before we embark on the detailed evolution of individual parameters, we will first 
consider the running of the quantities X1X2X3 and yiy2yz- They may be regarded as the 
analogs of phase space volumes and are measures of the overall magnitudes of the mixing 
parameters. We present the results in Fig.||, where we assume that, at high energy {t = 30, 
or E = lO^^GeV), the initial values are close to the "symmetric point" (xi = —yj = 1/6). 
We also assume that the Higgs Yukawa coupling constants are large (~ 0(1)), and not far 
from being degenerate. (The graphs are not sensitive to the precise initial values chosen in 
presenting the plots.) These choices, apart from esthetic reasons, are also motivated by the 
known neutrino parameters, which are usually considered to be a reflection of high energy 
behavior. It is seen that, with these initial conditions, Y\xi and Y\yj drop considerably as 
t decreases. However, the low energy physical values (n^i — 10~^, Y\yj — 10~^) require 
a drop ~ 10^^ from the symmetric point, and neither the SM or MSSM model can yield 
that. On the other hand, if we have some new model at high energies, it is possible that 
such a precipitous drop may be achieved. For this purpose we will simulate their effects by 
changing the values of b and c in Eqs.(3^, 3^, 3.21 , 3.22| ). This seems reasonable since, 
for the known examples, the basic structure of the RGE is the same, with the coefficients 
varying from one model to another. In addition, as we will see in Appendix ^ the two- 
loop RGE are also very similar and can be approximated, at least for a range of t values, 
by using effective b and c values in these equations. Another example can be found in 
extra dimension models. Although such models are non-renormalizable, it has been argued 
that the effect of introducing a cutoff amounts to scale dependent couplings which can be 
described by equations that are just like renormalization group equations [17, 18]. For 
the Yukawa matrices, in fact, the resulting equations are similar to Eqs.( |3.l| , p.2| ), where 
the coefficients b and c become proportional to the number of Kaluza-Klein states that 
contribute at a given energy scale. Thus, in our attempt to assess the behavior of the 
RGE, we will treat b and c as free parameters. As we will see, it turns out that, with 
our choice of initial values, the evolution does not vary much as b changed, but it is very 
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sensitive to the value of c. Thus, it is found that, for c ~ —7, Y\xi and Y\yj can indeed 
change dramaticahy in the neighborhood of t ~ 30. 

Turning to the individual parameters, we present in Figs.(^, ^ the evolution of (xj, yj) 
and the mass ratios. Using the same set of input initial values {xi ~ 1/6, yj ~ —1/6), we see 
that they behave similarly to products Y\ Xi and Y\ Vj^ albeit with slower rates of evolution. 
For the {xi,yj) parameters, with c = —7, we find that the evolution occurs rapidly and 
almost saturates within At ~ 5 — 10 from the initial t value, with xi — > 1 and all the other 
{xi,yj) 0. For c = — 3/2(SM) and c = l(MSSM), we do not find such convergence. The 
evolution of the mass ratios behaves similarly, although even with c = — 7, the hierarchy 
near t = does not approach the magnitude observed experimentally. Also, changing 
the b values has only minimal infiuence on these general results. Despite the similarity of 
Eqs.( |3.2i| ) and ( 3.22 ) to Eqs.( |3.8| ), ( |3.9| ) and ( 3.12| ), the numerical rates of evolution for 
the mass ratios differ considerably from those for the mixing parameters. A possible cause 
may be our choice of initial input values. However, so far our search for a region of input 
parameters which may result in a different behavior has not been successful. 

In summary, we have used numerical analysis to explore the properties of the RGE 
for the {xi,yj) parameters and the mass ratios. We assume that, at high energies, all of 
the (xj, yj) values and the coupling constants (/, h) are 0(1), and that (A/, Ah) are small. 
We also take the coefficients b and c as free parameters. We find a physically interesting 
scenario if c is negative and large. In this case, hierarchies in (x, y) develop rapidly at high 
energies. It is expected that, at intermediate or low energies, the model will be superseded 
by other models (such as the SM or MSSM). However, the hierarchical values are close to 
a fixed point of the RGE, so we will not see much evolution in these energy regions. Thus, 
the observed hierarchy is primarily accomplished at high energies, and will receive at most 
minor corrections at lower energies. At the same time, the hierarchies developed for the 
mass ratios are still not large enough to match those of the physical values. So the scenario 
of "infrared hierarchy" is only partially fulfilled. Further work is needed in order for one 
to arrive at a realistic model where the observed values of both (x, y) and mass ratios can 
be reached within this framework. 
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Figure 1: The evolutions of (a) xiX2Xs, (b) yiy2V?, and (c) ,P for SM(thin), MSSM(dashed) and 
the model(thick) with (6, c) = (3, —7), with the boundary conditions [{xi,X2) = (1/6 + e, 1/6 — e), 
(2/1,2/2) = (-l/6 + e,-l/6 + e), e = 0.01] and [fi = (1.6,1.8,2.0), hi = (0.5,0.7,0.9)] at t=30. 




Figure 2: The evolutions of (a) Xi and (b) for SM(thin), MSSM(dashed) and the model(thick) 

with (&, c) = (3, —7), with the boundary conditions [(.xi, X2) = (1/6 + e, 1/6 — e), {yi,y2) = (~l/6 + 
e, -1/6 + e), e = 0.01] and [fi = (1.6,1.8,2.0), hi = (0.5,0.7,0.9)] at t=30. 
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Figure 3: The evolutions of (a) /i/./, and (b) hjhj for SM(thin), MSSM(dashcd) and the 
model(thick) with (6, c) = (3,-7), with the boundary conditions [{xi,X2) = (1/6 + e, 1/6 — e), 
(2/1,2/2) = (-l/6 + e,-l/6 + e), e = 0.01] and [fi = (1.6,1.8,2.0), hi = (0.5,0.7,0.9)] at t=30. 
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5. Conclusion 



In this paper we studied the RGE of the quark mass matrices, using the recently proposed 
flavor mixing parametrization which is manifestly rephasing invariant. The resulting set 
of equations is not as complicated as the one using different parametrizations. Also, the 
dependence on the variables entering the RGE is quite transparent and simple. We also 
presented, for some simple functions of mass ratios, evolution equations which bear close 
resemblance to those of the {x, y) parameters. Since all of these quantities can be expressed 
as powers of A^, it is tempting to theorize that the hierarchies observed in mass ratios and 
in (x, y) may both originate from renormalization effects. In fact, when we substitute their 
known physical values into RGE, it is seen that one is close to a fixed point and there 
is little evolution for these parameters. Given the uncertainties in the initial values and 
the theoretical models, extrapolation from low energies is unlikely to yield an accurate 
picture of renormalization effects. Rather, to assess the general nature of the RGE, it 
seems more appropriate to start from a point with fast evolution, so that most changes 
are accomplished in its neighborhood, with minor corrections afterwards. To this end we 
assume that, at high energies, the quark masses are large but nearly degenerate while (x, y) 
values are not far from the "symmetric point", Xi = —yj = 1/6. It is found that, for a 
range of initial values, the parameters do evolve toward hierarchy. However, to reach the 
observed large hierarchy, one has to invoke the existence of new theories which give rise to 
large c values contained in Eqs.( |3.l| , f3.2| ). On the other hand, even with these parameter 
choices, the mass ratios do not match the large hierarchy in the observed values. Thus, we 
may have established an outline of an "infrared hierarchy" scenario. The details, however, 
are still lacking. Hopefully, our analysis can provide the impetus for further researches 
along this direction. 



Appendices 

A. Properties of RGE Matrices 



In Eqs. (|3.8| , |3.9| ), the RGE contain the matrices {Ai, Bi, A[, B[), which are listed explicitly 
m Table || We will now discuss some of their properties in detail. 

We note first that, since ^ ^fij = A/ijj = 0, the RGE are invariant if a constant is 
added to each column of [Ai, Bi, A[, B[). For instance. 



Ai > Ai + 



/ h h h 

h 62 63 I (A.l) 



\5i 62 53 

leaves Eq.( p.8| ) invariant. The symmetric pattern of the matrices listed in Table ^ is thus 



not inherent and can be changed by transformations as in Eg. ( A.l ). The constraints on the 
(x, y) parameters imply the existence of consistency relations among the {A, B) matrices. 
Indeed, it is found that X](^j ~ = ^i^i — B'^ = 0, which ensures Xi — Y^ yi) = 0. 



However, these relations are valid only if in Eq.( [A.l| ) all 5j = 0, j = 1, 2, 3. To account for 
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the arbitrary choice of 6j, we have the more general relations: 

(1 - P) Y^iA, - A'i) =0, (1 - P) Y,{Bi - Bl) = 0. 
Here, P is the cyclic permutation operator 

/o r 



(A.2) 



P 



1 I , 
1 



(A.3) 



so that (1-P) annihilates the 5 matrix in Eg. ([01), and (A/23, A/31, A/12) = (/I; /I; /i 
P). One also establishes the relations 

(1 - P) ^[(x, + Xj)Ak - {y^ + yj)^,] = Y^^XiXj - y,yj){l - P)w, (A.4) 

eye i>j 
(1 - P) J^[(x, + Xj)Bk - {y^ + yj)B'^] = ^^(x.x,- - y,y,)(l - P)w'^ . (A.5) 

eye i>j 

They validate the constraint ^[^(xjXj — yiyj)] = 0. When we compute ^(xiX2X3) and 

i>j 

^{yiyiyz), we find 



(1 - P) Y,{xiXjAk - y^yjAi) = 2j\l - P)w, 

eye 

(1 - P) Y^ix^XjBk - yiyjB'k) = 2j\l - P)uF 



(A.6) 
(A.7) 



eye 



If we put these two equations together, we obtain the evolution equation for J^, Eq.( 3.12| ), 
given in Sec.^. In the above equations, ^ denotes summation over cyclic permutation of 



eye 



{i,j,k) SO that, e.g., ^{xi + Xj)Ak = (xi +^2)^3 + (a;2 + a^s)^! + (a^s + a;i)^2, J^^iXjAk 

eye eye 

xiX2^3 + X2X3A1 + X3X1A2, etc. 

The matrices {Ai, Bi, A[, B'j) are also related to the w matrix. Let us write 



w = X + Y, 



where 





^ Xi X2 X3 ^ 




( Vi 


2/2 




X = 


X3 Xi X2 




2/2 


2/3 


yi 




X2 X3 Xi j 






yi 


y2 j 



(A.9) 



Then, e.g., the matrix A2 can be written as 





^ xi y2 xz\ 


^ X2 







^ yi X2 ys \ 


^ yz 





o\ 


A2 = 


x-i xi yi 





X2 




y2 ys X2 





y2 







\ ys X3 xi J 


^0 





X2 J 


\X2 yi y2 ) 


^0 





yi J 



(A.IO) 



So, to obtain A2, which is contained in ^X2, we begin by exchanging the X2 entries in X 
with corresponding ones in Y , and then multiply these matrices by diagonal ones composed 
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of the exchanged elements in X and Y . Similar rules apply to the construction of the other 
matrices. Notice also that the matrices are related by exchanging X2 ^ x^: 



Ai ^ Bi; A2 ^ B3; A3 ^ B2; 
A[ ^ B[- A'^ ^ B'^- A'^ ^ Bi,. 



(A.ll) 



B. Two-loop Renormailzation 



The one-loop RGE, Eqs.(3.1, |3.2D , must be amended when higher order effects are consid- 
ered. In particular, two-loop contributions can be included by writing 



(B.l) 
(B.2) 

where D^^^ are the one-loop contributions as given in Eqs.(|3.l|, |3.2| ), and D^^^^ are given by 



,(2) , J2): 



.(2) 



7(2) ^J2 



(2) 



.(2); 



D 



(2) 



d - {«?^ + b^Md + cfMu + dfMl + e'^'MdMu + f/>M^Md + gfUj, Md}. 



.(2) 



7(2) 



,(2) 



.(2); 



.(2) A#2 



(B.3) 

The coefficients a through g are given explicitly in Ref. [9] and will not be reproduced 
here. 

Following the same procedure as in the one-loop calculations, we find the RGE for the 
(x, y) parameters: 



-^(2)^;. = Af.Ai.H + Ah.Bi.F + Af.Ai.H + Ah.B-i.F, 



(B.4) 



-^(2)y. = Af.A'i.H + Ah.B'i.F + Af.A'i.H + Ah.Bi.F, 



where 



J' 

2el2)(/?/|) + + ff)\ /A/,„ H,, = [2e^^\hjh]) + ff\ht + /Ah,,. 

Note that ^ ^fij = A/ij^ = 0. Also, the structure of these equations is similar 
to that in Eqs.(3.8, 3.9). In addition, the matrices {Ai,Ai^Bi,B'j) enter these equations 
sandwiched between functions whose dependence on A/jj and A/ijj is the same as in 



(B.5) 



A/., = [cf + df{ff + /|)]A/,„ All,, = [c(2) + 42) (^2 + h])]Ah,,, (B.6) 



Eqs.(|3.8|, 3.9). Thus, the effect of two-loop contributions may be approximated by a change 
in the value c in Eqs. (|3.8| , ^ ), especially in the neighborhood of degenerate masses. 
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